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Abstract 



Last three years have seen new developments in the theory of last passage percolation, which 

has variety applications to random permutations, random growth and random vicious walks. It 

turns out that a few class of models have determinant formulas for the probability distribution, 

^ ' which can be analyzed asymptotically. One of the tools for the asymptotic analysis has been the 

l~/ ' Riemann-Hilbert method. In this paper, we survey the use of Riemann-Hilbert method in the last 

f^ ' passage percolation problems. 

o 

(^ ! 1 introduction 

-(— > 

s 

>'■ ifiiz) ^ e'^'+'"\ (1.1) 

j_j ■ Consider the following Riemann-Hilbert problem : Y{z) is the 2x2 matrix-valued function satisfying 

Y(z) is analytic in z G D,±, and is coninous in Q±, 

(1 Z~^Lp{z) , 
' '|, ^es, (1.2) 



Let S be the unit circle in C, oriented counterclockwise, and let il+ = {z e C : |z| < 1}, 51- = {z e 
\z\ > 1}. Set 



y(z) = (/ + 0(z-i))z''^'"-% asz- 



CX). 



Here CT3 = ( J _?i), and l + (z) (resp., Y-{z)) is the limit from the inside (resp., outside) of E. 

This Riemann-Hilbert problem appears in a class of problems which have many different interpreta- 
tions like longest increasing subsequence in a random permutation, last passage percolation, polynuclear 
growth model, and random vicious walks. In this paper, we survey a list of problems related to the 



above Riemann-Hilbert problems and discuss various aspects of the RHP (1.2) 



Before we state the problems, we first consider the solution to (1.2). From the work of Fokas, Its 



and Kitaev |20(], the RHP (O) is related to the orthogonal polynomials on the unit circle. Indeed (see 
e.g, Q), if we let 7rfc(z) be the fc*'' monic (leading coefficient=l) orthogonal polynomial with respect to 
the measure ip{z)dz/{2T:iz) on S, 

dz 

TTt^{z)z3^{z)—— = NkSjk 0<j<k, (1.3) 

|z| = l ^TTIZ 



the function 






is a solution to (1.2), where tt* (z) = z"7r„(l/z), and (C/)(z) is the Cauchy transform of / : 

Also it is a standard argument to show that the solution is unique. Here if can be a general function, 



not necessarily of the form (1.1). Thus for instance, 

TTkiz) = Yn{z; k), N-\ = -y2i(0; k). (1.6) 

In this paper, we discuss three different uses of this RHP formulation for orthogonal polynomials 
TTfe : (1) obtain asymptotic result using the Deift-Zhou steepest-descent method (2) obtain differential 
(or difference) equations (3) make a connection to so-called integrable operators. These three topics 
are discussed in Section ^, Section || and Section |^, respectively. The rest of the paper is organized as 
follows. In Section S, we state the last passage percolation problems whose solution has Toeplitz/Hankel 
determinant formulas, which can be expressed in terms of the solution to the RHP ( |l.2| ). These percola- 
tion problems have various different interpretations, and Section discusses some of these applications. 
There are a few other last passage percolation problems which can also be solved in terms of the RHP of 
the type (1.2), but with a different function tp. These problem are discussed in the final section ^. Most 



of the material in this paper have appeared somewhere. We indicate the references in each section. The 



difference equation. Theorem 5.1 is the only new result. 
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2 Last passage percolation 

Square case 

Consider a Poisson process of rate 1 in the plane. An up/right path is, by definition, a piece-wise linear 
curve with positive slope where defined. The 'length' of an up/right path is defined by the number of 



Poisson points on the path. Let L{t) be the length of the 'longest' up/right path starting from (0,0) 
ending at {t, t) (see Figure ^. See Section |^ below for other different interpretations. The basic formula 




(t,t) 



(0,0) 



Figure 1: Poisson points in (0,0) x {t,t) and the longest up/right path 



for us is the following, proved by p^ (see also [p8[) 
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^m<i)^-Eueuie){detMU))) 



(2.1) 



where (p is given in (1.1) and Z — e* is the normalization constant to make the left hand side 1 as 

As is well-known, using the Weyl's integration formula for the unitary group, the above expected 
value is equal to a Toeplitz determinant : 



Ec/gc/(£) det{(p{U)) = det{ipj^k)o<].k<e =: Di{t), 



(2.2) 



where ipj is the j Fourier coefficient of (p. Now we will relate this Toeplitz determinant with the RHP 



(1.2) 



As in Section]^, let 7rfc(z) be the fc*'' monic orthogonal polynomial with respect to ip{z)dz/{2'Kiz) on 



the unit circle and Nk be defined by (1.3). It is direct to check that the orthogonal polynomials have 
the determinant expression (see e.g., p^ ) 



TTkiz) :-- 



Dkiip) 



det 



1 ^0 
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•p-k+1 


<P-fc+2 • 


ipi 
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■ ^' J 



(2.3) 



Thus we have 



Ni, 



Dk+i 



(2.4) 



Note from (2.1) and (2.2) that 



Z = lim Dn, 



(2.5) 



which also can be seen from the strong Szego hniit theorem. Hence we have from (1.6) 



nm<i)^e~'"D,= lim ^ = []7V-i = [](_r2,(0;fc + l)). (2.6) 






Remark 2.1. In addition to the Toeplitz determinant formula, there are two other determinant formulas 
for¥{L{t) < £). Both of them are Fredholm determinants of operators, one acting on S (see /^), and 
the other acting on the discrete set {£,£+!,■■■} (see lljj , \2{\ , [Sj/j. The first Fredholm, determinant 
has the RHP expression which is algebraically equivalent to (1.2). 5*66 S'ection below for details. The 



second Fredholm determinant also has the RHP expression, but now the RHP has a jump condition on the 
discrete set. However, for the second formula, the kernel of the operator has integral representation, and 
hence the classical steepest- descent method is sufficient to obtain the asymptotics (US, W^J). Thus the 
lim,it of¥(L(t) < £) can be obtained from the second algebraic formula, which avoids the RHP analysis 
in Section^ below. But for some applications (see \28, U^) one also needs uniform, tail estimates of the 



'scaled' random variable {L{t) ~ 2t)/t^'^ in t (see ( 4.17 ); ( 4.18| ) below). This problem is more difficult 



than the convergence in distribution, and the lower bound (4.18) has been obtained only from RHP 
analysis so far J^. The identities between the Toeplitz determinants and the Fredholm determinants are 
discussed in ^, |7|, |7|, ||, [7|, \l^. 

Triangle case 

Suppose now that we take a (2-dinicnsional) Poisson process of rate 1 in the half plane y < x, and 
take a (l-dimensional) Poisson process of rate a > on the line y — x (see Figure ||). Let Ls{t; a) be 

(t,t) 



ro.m 

Figure 2: Triangle case 

again the length of the longest up/right path from (0, 0) to (i, t). The algebraic formula for this case is 
obtained in pq] for a = and in H for general a > : 

P(L,(i; a)<£) = ^ ^ueow det((l + at/)e*^), (2.7) 

where Z = e"*+*'/2. 

There are two components of 0{i) depending on the sign of the determinant of the matrix. Using 
the Weyl's integration formula, the expected value over each component is expressed in terms of one of 



the Hankel determinants of the form (see e.g., Theorem 2.2 of 0) Hi = dct{hjj^k)t-><i where 



hj^ j a;J7i(2:)(l-a;)±2(i + 2:)±2rf2;, (2.8) 



ith 



1 



h{x) = {\+a^ + 2ax)e^*'' ^{l + az){l + a/z)ip{z), x= -{z + z^^). (2.9) 



As an analogue of (2.3), these Hankel determinants are related to the orthogonal polynomials on the 
interval (—1, 1). But the orthogonal polynomials on the unit circle with weight (1 + az){l + a/ z)ip{z) 
and the orthogonal polynomials on the interval (—1, 1) with weight h{x){\ — a;)* 2 (l + x)* 2 are related 



in a simple way (see 43 ). Thus the ration of Hankel determinants Hi/Hi^i can be expressed (see 
theorem 2.3 of Q) in terms of orthogonal polynomials (norms and the values at z = 0) on the unit 
circle E with respect to the measure 

(l + az)(l + a/z)^(z)-^. (2.10) 

ZTTJZ 

Moreover, after some algebraic work (see Section 3 of H), the simple factor (l+az)(l+Q;/z) in the above 
measure can be removed and Hg/ Hi^i can be expressed purely in terms of orthogonal polynomials ttj, 
for the measure ip{z)dz/{2'Kiz). The result is ^ 

nL,{t;a) < 2^+ 1) = e-"*i{(7r*^(-a) + aTTa^-a))/// + (7r*^(-a) - aTT2i{~a))H-} (2.11) 



where 



H^ = n ^2fc+i(l T 7r2fc+i(0)). (2.12) 



k=i 



There is a similar formula for ¥{Lt{t; a) < 2£). Thus again, the probability distribution for Ls{t; a) can 



be expressed in terms of the solution of the RHP (1.2) 



Remark 2.2. In /0, ^a, W, the authors considered five types of symmetry of the Poisson model. Re- 
flection symmetry about the anti-diagonal, reflection symmetry about both diagonal and anti- diagonal, 
and rotation symmetry about the center are considered in addition to the square case and the triangle 
case. The distribution of the longest up/right path in each case has the determinantal expression (either 
Toeplitz or Hankel determinants with a simple change of the weight), and can be expressed in terms of 



the RHP (1.2 



Remark 2.3. As in square case (see Remark 2.1), there is also a different algebraic formula for the 
triangle case /pT/. But the result is so-called Fredholm Pfaffian of an operator acting on a discrete set. 
The kernel of the operator is rather involved, and so far there has been no work for the asymptotic 
analysis from this Fredholm Pfaffian formula. 



External sources 

Suppose that in the square model above, we have additional (1-diniensional) Poisson processes of rate 
a+ and a_ on the lines y — and a; = 0, respectively. We assume that the corner (0,0) has no point 
(see Figure 0). Let Le(t;a+,a_) be the length of the longest up/right path from (0,0) to {t,t). This 



(0,0) 



-^ (t,t) 



Figure 3: External sources case 

problem arises from a polynuclear growth model with random initial data 
The distribution in this case is given in H] : 



P(Le(i; a+, a-) = -{D'^ - a+a-D[_{) 



where 



(2.13) 



D', := Ec/ec/(£) det((l + a+U){l + a_C/-i)^(C/)), 



and z = e"+*+"-*+* 



By a similar argument as in (2.10), D'^ can be written as (see Theorem 3.2 of H) 



, _ 7r^(-a+)7r*(-a_) - a+a-TT(,{-a+)TTt{-a^) 



D'e = 



Dp, 



a^a- 7^ 1, 



(2.14) 



(2.15) 



1 — a+a- 

where Dp = Di{t) is given in ( ^.2| ). For a+a_ = 1, we use the I'Hopital's rule in the above formula. 
Thus again the distribution for Le(i; a+, a-) can be expressed in terms of the solution Y of the RHP 



(1.2) 



Remark 2.4- In the triangular case above, we may put an additional (1- dimensional) Poisson process of 
rate a+ on the line y = 0. Then there still is a determinant formula for the distribution of the longest 
up/right path, which again can be expressed in terms of the RHP (1.2). 



3 Applications 



In this section, we discuss various different interpretations of the last passage percolation problems in 
Section 0. 



Random permutations 

Given a permutation tt, a subsequence {7r(ii), • • • , 7r(ife)} such that ii < ■ ■ ■ < i^ and 7r(ii) < • • • Tr{ik) 
is called 'increasing'. The length of such an increasing subsequence is defined to be k. Let Ln{tt) be 
the length of the longest increasing subsequence of a permutation n e Sn- Now take tt randomly from 
Sn, and we ask the distribution of L^. This problem is known as Ulani's problem since early 1960's 
(see e.g., @, |). 

If we take N itself as a Poisson random variable of rate t^, the distribution of Ljv is same as the 
distribution of L(t) above : 

oo —t^ (±2\N 

nm < ^) = E ' jv! ^^^^ - ^)- ^^-^^ 

Af=0 

In other words, f'{L{t) < £) is the exponential generating function of P(LAr < £). This follows from the 
fact that a configuration of N points in a square can be regarded as a permutation tt G Sn '■ consider 
the relative a;-orders, then relative y-orders. The longest up/right path in a points configuration is 
precisely the longest increasing subsequence of the corresponding permutation. Thus the asymptotics 
of Ln as N ^ oo can be obtained from the asymptotics of L(t) as t — » cxd (see pq). 



The limiting distribution of L(t) (see (4.16) below), and thus the limiting distribution of Ln, is 



first obtained in [^ by using the RHP (L2). Later the same result was obtained by ||3^, [l^, ^ using 



different methods which avoids the RHP analysis. The uniform tail estimates (see ([4.17), (|4.1§|)) and the 



convergence of moments of the scaled random variable were obtained in [^ . As mentioned in Remark 



2.1, the (lower) tail estimate obtained in [^ was later used in other applications |g8[, |4l[ |. 

Now we consider the corresponding permutation model for the triangle case. Given random points in 
the triangle model, take mirror image of the points in y < x about the line y = x. Hence we have points 
in the whole plane which is symmetric about y ^ x. It is each to see that in terms of permutation, the 
symmetry condition implies that tt'^ — 1. Also the points on the line y = x corresponds to the fixed points 
7r(i) = i. Thus it is natural to consider the set Sm.n — {tt & S2m+n : tt^ = 1, \{x : tt{x) = x}\ = n}. We 
equip this set with the uniform probability measure and consider the longest increasing subsequence 
of a random tt e Sm,n- As before, the asymptotics of Ls{t]a) implies the asymptotics of the longest 
increasing subsequence in this set. See [g|, |l^, ^ for details. 

Plancherel measure 

Let Yn be the set of partitions (Young diagrams) A = (Ai > A2 > • • • > Xn), Aj G NU {0}, ^ ■ Xj = N. 
Given a partition A, let d\ be the number of standard Young tableaux of shape A (see e.g. [^ for the 
definitions). It is a basic fact in the representation of symmetric group that X^ahw ^a ~ l"^^! ~ -^■• 
Define the Plancherel measure on Yn by 

P(A) = ^, XeYN. (3.2) 

The famous Robinson-Schensted (RS) correspondence (see e.g. |3^, |Q) implies that the first part 
Ai from the Plancherel measure has the same distribution as the longest increasing subsequence of 



random permutation tt G Sn from the uniform measure. Thus from the resuh of the fongest increasing 
subsequence, the asymptotics of the first row Ai can be obtained. The convergence in distribution and 
the convergence of moments for the second row A2 was obtained in M. For the general rows Aj, j > 3, 
the convergence in distribution is obtained in |3^, |l3|, ^ , and the convergence of moments are obtained 
recently in B. 

Let /9 > 0. As a generalization, consider the /?-Plancherel measure on Y^ defined by 



P(A) = ' p , A e Yn. (3.3) 



One can regards /3 as the inverse temperature. The above Plancherel measure corresponds to the case 
when (3 = 2. The 1-Plancherel measure also has a simple permutation interpretation. Consider the 
set of involutions Sn '■= {tt G Sn ■ tt^ = 1}. Equip Sn with the uniform probability measure and 
consider the length of the longest increasing subsequence of a random vr G Sn ■ The RS correspondence 
implies that the first row Ai of A G Yn under the 1-Plancherel measure has the same distribution as the 
longest increasing subsequence of a random involution. In ||l^ (see also |^), the limiting distribution of 
the first row of the 1-Plancherel measure was obtained. It would be interesting to obtain the limiting 
distributions for the general (3. 

Polynuclear growth model 

Consider a flat 1-dimensional substrate. There are random nucleation events happening on the substrate, 
which can be thought as a 2-dimensional Poisson process in space-time coordinates. An island of height 
1 with width is made at each nucleation position, and the island grows laterally with speed 1. When 
two islands meet, they just form one island and keep growing at the edges. Further nucleation events 
may happen on top of islands. This model is called Polynuclear growth (PNG) model. An interest is 
the height h{x, t) as time t tends to infinity. 

It is an observation by Prahofer and Spohn that for certain choices of initial data, PNG models are 
in bijection to the last passage percolation models in Section g. As the first case, suppose that initially 
there is an island of width at the position x = 0. It grows with speed 1, and we assume that further 
nucleation events occur only on top of the base island and on top of islands on the base island. Then the 
height h{0,t) is equal, in the sense of distribution, to L(i/-\/2) of the square case (see [£5|). Secondly, 
we consider the half infinite case : x G [0,oo). When an island hits the edge a; = 0, it stops growing 
at that side. As before there is the base island at the position a; = at time t = 0. In addition to 
the nucleation events of rate 1 in (0, 00), and we assume that there are extra nucleation events at the 
position a; = of rate a. Then the height h{0,t) is equal to L(i/\/2; a) of the triangle case (see |3^). 
If one considers flat case, i.e. no base island, the PNG model is in bijection to yet a different symmetry 



model of the Poisson process mentioned in Remark 2.2 above (see |35|). Finally the case of external 



sources in Section has the PNG model which has random initial data. See |R4, BSl for more details. 



Random turn vicious walks 

At time t = 0, each site of the semi-infinite lattice {• • • , —2, — f , 0} is occupied by a walker. A walker 
is called 'right-movable' if its right-neighboring site is vacant. At each discrete time, we pick a random 
walker among all right-movable walkers (there are finitely many such walkers) , and move it to its right- 
neighboring site. Hence at each time, one and only one walker moves to its right. We let the walkers 
walk by the above rule for N time steps. In the next N time steps, we let walkers move to their left 
according to the similar rule. We further impose the condition that after total 2N time steps, all the 
walkers return to their original positions. The interest is the position xi{N) of the right-most walker 
at time N. It is shown in |M (see also p4|) that the distribution of a;i(A^) is same as the distribution 
of Ln, the length of the longest increasing subsequence of a random permutation. 

If we just let the walkers walk up to time N without assuming that they all return to their original 
positions, xi (N) is, in the sense of distribution, equal to the length of the longest increasing subsequence 
of a random involution (tt^ — 1). See S for the asymptotic results of this case. 

Queueing theory and Totally asymmetric simple exclusion process 

Lattice versions of the last passage percolation problems are known to be related to the queueing theory 
and also so-called totally asymmetric simple exclusion processes (TASEP) (see e.g. |^). See Section ^ 
below for lattice last passage percolation problems. 

4 Asymptotics 



In this section, we sketch asymptotic analysis of the RHP (1.2) and discuss the asymptotic results for 
the random variables L{t), Ls{t;a) and Le(i; «+, «_). 



The asymptotics of the RHP (1.2) as k,t — > oo in all different regimes was first considered in 



and later in a little more detail in ||l^ . A special interest is the case when k and t are related as 



k ~ 2i/3fc2/3- y^-^) 



For X in a compact set, the RHP (1.2) is localized around z = — lasfc— >oo, and the local RHP is, 
after scaling, the RHP for the Painlcve II solution. In the below, we will discuss, heuristically, how the 
Painlcvc II solution arises. This heuristic argument is taken from Section 6 of Iq]. 

First, we show why the point z = —1 plays an important role in the asymptotic analysis. Since the 



RHP (1.2) is not normalized at oo, we algebraically transform it to a normalized problem. Set 

and set 

/q -1\ / z~'^e*^^\ 



Then m(z) solves the new Riemann-Hilbert problem 

m{z) is analytic in z G fl±, and is continous in fl±, 

m+(z) = m_(z)y(z), z e E, (4.4) 

Tn{z) — I + 0{z~^), z ^ oo, 

where 

V(.)=f .V --'^f(^)-^V ^W-e*^--). (4.5) 




Now the oscillatory factors in the jump matrix are V12 and V21 = — I/V12. It is direct to check that the 
critical points of 

logF2i(z) = fc(^(z-z-i) + logz) (4.6) 

are 



When 2t < k, there are two real critical points. These two points collapses to one point z = — 1 when 
2t = k. And when 2i > fc, there are two complex critical points. Thus one can imagine that the 
nature of asymptotics changes at 2t — k. This is indeed the case. It is well-known that in the (so- 
called Plancherel-Rotach type) asymptotics of the orthogonal polynomials, the equilibrium measure of 



an associated variational problem plays a key role (see e.g. |17| [l^]). For the case at hand, the support 
of the equilibrium measure is the full circle when 2t < k, and it is a part of the circle when 2t > k. And 
the point z = — 1 is the place where the 'gap' of the support starts to open up at 2t — k. 



Thus it is natural to analyze the problem near z = — 1 and in the regime 2t ^ k. The scaling (4.1) 
is chosen so that the local RHP around z = — 1 becomes non-trivial. Writing z = — 1 + ,s for z near —1, 
we have 

logy,2(.) - m = -^,{k^''s) - ^^ik'^'s)^ + ^^(1 - ^ttIst,) + • • • . (4.8) 

If we take w — 2^^^^k^^^s^ then we have 

log 1^12 (z) - 7ri = ~2xw + - w^ H , (4.9) 

and hence the jump matrix near z — —1 becomes as fc — > 00, 

for w E iM.. But after rotation by 7r/2 and removing (—1)'" by a simple conjugation, this is precisely the 
RHP for Painleve II equation with the choice of parameters p = —q = 1, r = (see e.g. [E9|, 0). Of 

10 



course, this is only a heuristics, and one needs to justify the convergence of the original RHP (1.2) to 
the above local RHP. As mentioned above, the justification proceeds differently for 2t < k and 2t > k, 
due to the change of the support of the equilibrium measure. See Q, [Q for details. 

The asymptotic results we obtain are, for example, the following [Q, |l^ : There are numerical 
constants C, c, xq > such that for large k and t, 



(i). if a; > xq, 



(ii). if —xq < X < Xq, 



-Y2ii0;k)-1 



Yu{0;k) 



< 



C 



(iii). if X < — xo. 



2i/3 
-Y2i{0;k)-l-^^v{x) 



k(2t/k-log{2t/k)-l)Y^^fQ. fc) _ 1 



yn(0;fc) + (-l)'=^«(x) 



< 



_C_ 



(-1)' 



2t 



2t-k 



Yii{0;k) 



< 



C 
2t~k' 



(4.11) 



(4.12) 



(4.13) 



Here w(x), u{x) in case (ii) are given by the Painleve II solution 

u" — 2ii^ + xu, 

i_ 

2^x 



^(^) Sy^ryre ^^^"'' '> a; ^ +oo. 



and 



v{x) = / {u{s)fds. 



(4.14) 



(4.15) 



We also obtain the asymptotics of Y{z) for general points z e C |10| . 

Now using ( 4.1lD -(4.13), (2.6) and (1.6) yield the convergence in distribution in the square case 
: for fixed a: e M, 



y^U^^^ < x ) = exp( - / (y - x){u{y)fdy ] =: Fgue{x). 



(4.16) 



The function Fgue{x) in ( 4.16 ) is called the GUE Tracy- Widom distribution [^, which was first 
obtained as the limiting distribution for the fluctuation of the largest eigenvalue of a random matrix 
from the Gaussian unitary ensemble in the random matrix theory (see e.g. ||30|1 for an introduction 
to random matrix theory). Also the uniform tail estimates can be obtained M : there are constants 
C, c,xq > such that for all i > 0, 

'L{t)-2t 



and 



il/3 



L{t) - 2t 

il73 



<x\ <Ce' 



<x\ <Ce~ 



X > Xq, 



X < —Xq- 



(4.17) 



(4.18) 
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The upper estimate (4.17) also follows from a large deviation result of BO]. But the lower estimate (4.18) 
is obtained only from the RHP analysis so far. These estimates imply the convergence of moments of 
the scaled random variable {L{t) — 2t)/t^'^ [Q. These estimates have also been crucial for the analysis of 



the transversal fluctuation of the longest up/right path 1 28 and of the perturbation of the equilibrium 
measure for a related dynamical system, called stick process |^ . 

The asymptotic analysis of the RHP (1.2) also yields similar results for the triangle case Ls{t] a) pO|| 



and the external sources case Le(i;a+,a_) 0. The results show more interesting feature. Depending 
on the values a, a+,a_, the scaled random variables converge to different distributions. For example, 
{Ls{t\a) — 2t)/t^/^ converges in distribution to the so-called GSE TW-fmiction [|5|, 

FasEix) := i (ei ■''" ^(y)dy + ^-h\T «fe)'^A (FauEix))^''' (4.19) 

for fixed < a < 1, and to the GOE TW-function g|], 

Fgoe(x) := eh ^^ ^^■y'>''\FGUE(x)fl^ (4.20) 

for a = 1. When a > 1 and fixed, a differently scaled random variable converges to the Gaussian 
distribution. Moreover if we scale a ~ \ — -|%- with fixed w g M and take i ^ oo, the above scaled 



random variable converges to yet another distribution function for each fixed w G M. This new one- 
parameter family of distributions interpolates Fqse and Fgoe as w — oo and w — Q. Similar feature 
appears for the analysis of Le{t; a+, a_) where two-parameter family of new distributions are obtained 

0- 

5 Difference equation 



In this section, we show that the RHP (1^) yields difference equations for the entries of Y{0). It is 



direct to check from (1.4) that yii(O) = — y22(0). Also since the determinant of the jump matrix for 



(1.2) is 1, we have detr(O) = 1. 



Theorem 5.1. Let 



r(0;.)= -''"^> ''"^> . (5.1) 

I a{k) bik) 



Then b satisfy the discrete Painleve II equation 



and a and d satisfy 



h{k) + {b{k - 1) + 6(fc + 1)) (1 - b{kf) = 0, (5.2) 



a{k) = (l-6(fc)2)a(fc + l), (5.3) 

dik) = (l-6(fc)2)d(fc-l). (5.4) 
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Proof. Consider (4.4). Note that the jump matrix has determinant 1. Thus det(r7i-|_) = det(?7i_), and 
hence det(?7i) is an entire function. Also det(?7i) ^ 1 as z -^ oo. Thus by LiouviUe's theorem, we have 
det(TO(z)) = 1, and in particular, (m{z))~^ is analytic in n±, and is continuous in fl±. 
Set 



m{z) = m{z: k) = I -\ ^-^ + 0(z"^), z~^oo, 



and set m(0; k) — B(k) and 

m{z) = m{z; k) = B{k){I + m\{k)z + 0{z^)), 
Set 

*{z) = *(z;fc) :=m(z) I ^ ,° ), Ji := ( ^ ^ 



z^O. 



Jn:= 



'1 0^ 
.0 0, 



zJi--^Jl 



Z —> OO, 



Note that by (|5J) and (^), 

*(z;/c) = S(fc)(/ + m°(A:)z + 0(z2))e-^-''i+'=^i°«^''i-*-^i+*^^\ z ^ 0. 
Now the jump matrix for ^(z) is a constant matrix : 

' l\ 



vl/+(z) = *_(z) 



-1 1, 



ze s. 



(5.5) 
(5.6) 

(5.7) 

(5.8) 
(5.9) 

(5.10) 



Equations : By differentiating (5.10), ^^ satisfies the same jump condition. Hence §7^"^ has 
no jump cross S. Also from (5/7), §f-^~^ has a double pole at z = 0. Therefore using (^.8| ), (x£), we 
have 






(5.11) 



9z z z" 

for some constant matrices Ai,A2 which depend on k. On the other hand, ^(z; k + l)\E'(z; fc)^-'^ again 



has no jump cross E, and it is now entire. Thus using ( |5.8| ), ( |5.9[) , we have 

*(z; fc + l)*(z; fc)~i = zJi + Xo(fc) =: g(z; fc), 

for some constant matrix Xo{k). 

Now we take 'cross differentiation' of (5.11) and (5.12|) in two different ways : 



(5.12) 



^^{z; fc + 1) = ^Q{z; k)^(z; k) = ^^^^{z; k) + Q{z- fc)^5l!^ 
az az oz oz 



\^^^^ + Q{z;k)P{z-k)Y{z;k). 



(5.13) 



and 



9*(z;fc) 5\E'(z;fc) 



S.2 



9z 



P{z; k + l)*(z; /c + 1) = P{z; k + l)(5(z; /i:)^(z; fc). 



(5.14) 



feH^fc+l 
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Thus we obtain an equation 



dQ{z;k) 
dz 



+ Q{z; k)P{z- k) = P{z- k + l)Q(z; k). 



(5.15) 



By plugging in the formulas of P and Q and comparing the coefficients in z, we obtain the relations for 



Ji + Ji^i(fc) - Ai(fc + 1) Ji + t[Xo, Ji] = 0, 

JiA2{k) - A2{k + 1) + Xo(fc)Ai(fc) - Ai{k + 1)Xo(.t) = 0, 

Xo{k)A2{k) - A2(fc + l)Xo(fc) = 0. 



(5.16) 

(5.17) 
(5.18) 



Constant matrices Ai, A2,Xq : Now we express Ai, A2, Xq in terms of m. We plug in (5.8) and 



(5.9) into (5.11) and (5.12). This determines the constant matrices Ai, A2,Xq. Using IpM for (5.11), 

0(z-'^) : Ai{k) = i[mf (fc), Ji] + kJi. 



Using (5.9) for (5.11) 



Using (5.8) for (5.12), 



O(z-i) : Ai(fc) = B{k){t[m'l{k), Ji] + kJi}B{k)-\ 
0{z-^) : A2{k) = tB{k)JiB{k)-\ 

0(1) : Xo{k) = Jo + m'^{k + 1) Ji - Jimf{k). 



Using dJ) for (|5.12D , 



0(z) : Ji = S(fc + 1){ Ji + m;(fc + 1) Jo - Jom1{k)}B{k)-\ 
0(1) ■.XQ{k)^B{k + l)JoB{ky\ 



(5.19) 

(5.20) 
(5.21) 

(5.22) 

(5.23) 
(5.24) 



(5.25) 



(5.26) 



Symmetry : Note that the jump matrix v{z) for m satisfies 

a3v{z~^yc73 = v{z), as-- I 

Thus it is standard to show that m satisfies the symmetry condition 

m{z) - a3B*(m(z"i)*)-V3. 

By taking z -^ 0, and considering terms of order 0(1) and 0{z), we obtain the symmetry relations 

B - fXaSVa, (5.27) 

m° = -asim^fas- (5.28) 

Since det B — 1 and tr{m1) = 0, we can set 

^ a(fc) 6(fc)\ 



B{k) 



-b{k) d{k) 



a{k)d{k) + b{ky ^ 1, 



(5.29) 
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and set 



^?(fc) := 



Mk) -A(fc), 



mf{k) := 



-A(fc) i^{k) 
fi{k) \{k)^ 



(5.30) 



Equations from ( 5.16| )-( 5.1S ) : Now we re-write the equations ( ^.16 )-( pTl8 ) in terms a,b,d and 
X,fj,,:y. First, if we use (5.19|) and (5.24), use (5.29|) and (5.3C), the equation (^.16) yields 



,j{k) = b{k + l)dik), 
j/(fc + l) = -a{k + l)b{k). 



(5.31) 

(5.32) 



Second, if we use ( 5.2C ), ( 5.21 ) and ( 5.29 ), ( 5.30| ), the equation ( 5.17 ) yields the same equations ( 5.31 ) 



and (5.32). Finally, the equation (5.18) is trivial if we use (5.21) and (5.24). 

Equations from ( 5.19 )-( 5.23| ) : We also have the condition that the two different formulas ( 5.19 ) 
and (5.2C) for Ai{k) are equal. Similarly ( 5.22| ) and ( 5.24| ) for Xo{k) are equal. Also there is an 
additional condition ( ^.23 ). From ( 5.19 ) and ( 5.2(]| ), we have 



(i[mf (fc), Ji] + kJi)B{k) = B(fc)(i[mf (fc), Ji] + fcJi). 
Using (HH), i ^M ), this yields 

t{n{k)a{k) - v{k)d{k)) + kb{k) = 0. 



If we use equations ( 5.31 ) and ( |5.32| ), and the relation ad + b"^ — \, this equation becomes 

kb{k) + t{b{k - 1) + b{k + 1)) (1 - b{kf)) = 0. 



From ( |5.22[ ) and (|5.24|) , we have 

(Jo + mf (fc + 1) Ji - Jim'^{k))B{k) = B{k + 1) Jq, 

which yields 

a{k+l)-a{k) + v{k + l)b{k) ^ 0, 

b{k) + v{k + \)d{k) = 0, 

b{k + 1) - {\{k + 1) - \{k))b{k) ^ ^i{k)a{k) = 0, 

(A(fc + l)-A(fc))d(fc)-/i(fc)6(fc) = 0. 

Thus using ( |53^ ), ( ^37| ) yields 

a(fc) = {l-b{kf)a{k + l). 

On the other hand, the equation ( 5.231 ) implies 

JiB{k) = B{k + l)(Ji + m\{k + 1) Jo - Jom\{k)), 



(5.33) 

(5.34) 

(5.35) 

(5.36) 

(5.37) 
(5.38) 
(5.39) 
(5.40) 

(5.41) 

(5.42) 
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which yields 



(A(/c + l)-A(fc))a(fc + l)-i/(fc + l)5(fc + l) = 0, 

b{k + 1) - n{k)a{k + I) = 0, 

b{k) - (Xik + 1) - X{k))b{k + I) - iy{k + l)dik + I) = 0, 

dik + 1) - a{k) + ^l{k + l)b{k) = 0. 



Thus using (5.31), (5.46) yields 



d{k) = {l~b{kf)dik-l). 



(5.43) 
(5.44) 
(5.45) 
(5.46) 

(5.47) 
D 



Remark 5.2. The RHP (1.2) has also the parameter t. In addition to the above equation (5.2)-(5.4) 
obtained from the z-derivative and 'n-derivative' of Y , we can also obtain two more equations from 
z-derivative and t-derivative of Y , and from the t-derivative and n-derivative of Y . In this paper, we 
would not consider those equations. See Section 3 of j^ and |^q , [(% |i|/ for those other equations. 



6 Integrable operators 



In this section, we related the RHP (L2) to a so-called integrable operator on L^(S), and also discuss 
an identity between Toeplitz determinants and Fredholm determinants on L^(E). 

It is shown in M that the RHP ( |l.2| ) for orthogonal polynomials is also the RHP for an integrable 
operator on E after a simple algebraic transformation. For this purpose, the particular choice (p{z) = 
^t{z+z ) jj^ (1.2) is not important. In the below, we just assume that f(z) is a continuous function 



on S which has the factorization ip{z) = ip^{z)(p^(z) where (p-^-{z),lp-^-(z)^^ are analytic in \z\ < 1, 
continuous in \z\ < 1, (p+{0) = 1, and (p-{z),ip-{z)^^ are analytic in \z\ > 1, continuous in \z\ > 1, 
liniz^oo </'-(^) = 1- We also assume that 1 is not in the spectrum of the operator K defined in (6.7) 
below. 



Define (cf. 




Yiz) 



i^+(z) 

, 2^+{z)-\ 



\z\<l, 



and set 



M{z) :-- 







Y{z) 





'z''ip^{zy 







\z\>l. 



(6.1) 



(6.2) 



Then M(z) solves the new Riemann-Hilbert problem 

M{z) is analytic in z G rt±, and is continous in n±, 

M+{z)^M-{z)V{z), zeS, 



(6.3) 



Af(z) = / + 0(z-i), 



CXJ, 
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where 

V{z) = (, ^^^ ^ -2^-V(^)-^\ ^^^^ ,^ ^+{z)^-{zr\ (6.4) 

Note that the jump matrix V is of the form 

Viz)=I-27Tif{z)g{zf, (6.5) 

where 

.f{z) = {fij2f = (^-^-i^(z))^, g{z) = {g,,g,f ^ -L(^fe, 2^(s)-i)^. (6.6) 

This is precisely the RHP canonicaUy associated to the so-caUed integrable operator K — Kk on L^(S) 
(see [^, |l^) whose kernel is defined by 

z — w 27n[z ^ w) Jy: 



It is shown in p3| (see also |16[) that given an integrable operator K of the form (6.7) (with general 
column vectors /, g such that f{z)'^g{z) = 0), the resolvent kernel of K , if exists, is again an integrable 
operator 

1 .^-),,„0 = a£)ISf) ,6.8, 



\ — K J z — w 

and moreover, the column vector F, G are expressed in terms of the solution M to the (normalized) 



RHP with the jump matrix (3.5) 



i^(^)=(Y^/)(^) = M^{z)f{z), (6.9) 

G(z)=(^^.g)(z) = (Af+(z)^)-\9(z). (6.10) 

From the structure of the integrable operators, one can also show (see Lemma 4 of [p|) that 

Thus, using the definition of Af , ( [Tq ) and ( ^.4| ), we obtain 

^fc-i(y^) ^ ^ det(l-jrfc„i) ^ 

i?fc((^) det(l-i^fe) ■ ^ ■ ^ 

Moreover, when [1, oo) has no intersection to the spectrum of Kk for all k > ko for some fco, it is shown 
in Lemma 5 of [^ that 

lim 2-P det(l - Kp) = 1. (6.13) 

p — >oo 
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Thus for this case, by taking infinite product in k of ( |6.12 ), we have 



where 



2-"det(l-X„), (6.14) 



£'ooM= lim i:'p(93) =exp(Vj(log(^)j(log(/7)_j) (6.15) 

p — >-oo ' ^ 

3 = 1 

is given by the strong Szego theorem. 

It would be of interest to find if the other entries M12, M21 have expressions involving K, analogous to 
(3.11). Such formulas, if exist, would yield similar formulas to (6.12), ( |6.14 ) for the Hankel determinants 



(2.12) 



We note that as mentioned in Remark 2.1, there is yet another identity between the Toeplitz de- 
terminant and the Fredholm determinant of an integrable operator. But in this case, the integrable 
operator acts on a discrete set {k, fc + 1, • • • }. Hence there is another RHP with jump conditions on the 
discrete set. It is not clear yet if for example, there is an algebraic transformations between these two 
RHP's, one with jump on E and the other with jump on the discrete set. 

7 Other models 

There are other last passage percolation models which have the similar Toeplitz/Hankel determinant 
formulas of the form (2.1), (2.7) and (p. 13). The only difference now is the symbol ip which depends 



on the model. In the below, we will describe the models and the corresponding symbols. Here we only 
consider the square case and the triangle case. 

Since the RHP formulation for the Toeplitz/Hankel determinants in Section || is general, for all the 
models below there are associated RHP of the form ( |l.2| ) with new function (p{z). Thus we can follows 
the procedures in Section |j. Section ^ and Section |6|. The determinant identity in Section |g is general 
and it includes all the cases below. The difference equations for the cases below can be worked out as 
in Section ^. It will appear in some other place. See Q for some results for the difference/differential 
equations. The asymptotic analysis of Section |j can be worked out in all cases below. However, as 
mentioned several places above, for the square case, the convergence in distribution result of the form 



( 4.16 ) can be obtained from a Fredholm determinant formula on a discrete set, which does not require 
the RHP analysis. In the literature, many such results are obtained from the Fredholm determinant 
analysis. But for the triangle cases and also for the tail estimates, RHP analysis have been used (see 
P, 0). For the survey of asymptotic results of the cases below, see e.g. p9|. 

7.1 Square case 

Consider a lattice version of the Poisson process. Let X{i, j), i, j € N, be a planar array of independent 
random variables. We consider a weakly-up/weakly-right path, which a sequence {(*fe, j*:)}fe^i such that 
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ii < ^2 < • • • < *; and ji < J2 < • • • < ji- Let (1, 1) /" (M, A^) denote the set of all weakly-up/weakly- 
right paths from (1, 1) to {M,N). Now define 

L(M,N):= max V X(i,j). (7.1) 

(ij)e7r 

If one thinks of the random variable X{i,j) as the passage time at the site {i,j), the sum J^a ,)g7r -^ih j) 
is the total passage time to travel from (1, 1) to {M, N) along the particular path tt. Thus L(M, N) can 
be regarded as the last passage percolation time from (1, 1) to {M, N). As in the Poisson case of Section 
0, this problem arises in many different fields. See for example, |27, ^ ^ for various applications. 



In the definition of L{M, N) in (7.1), we could take different admissible up/right paths. For example, 
we could take weakly-up/strictly-right paths (sequences {{ik,jk)}k=i such that ii < • • • < i;, ji < • • • < 
ji) or strictly-up/strictly-right paths (sequence {{ik,jk)}k=i such that ii < • • • < i;, ji < • • • < ji). Also 
we can take continuum model instead of lattice model. In the below, we list the models which have the 
determinant formula for the distribution of L{M, N). In each case, we have 

P(i(Af , N)<£) = -^ ^ueuw det(^(f/)). (7.2) 

'^ALN 

with a function ip{z) for \z\ — 1 which depends on the model. The constant Zm.n = lim^-too De(ip) is 
a finite number in each case. 

We use the notation g{q) for a geometric random variable with parameter q S (0, 1) : F{g{q) = 
k) = {1 — q)q^ , A: = 0, 1, 2, • • • . We understand that g{Q) means the identically random variable. The 
notation b{q) is used for the Bernoulfi random variable : V{b{q) = 0) = j^, P(6(g) = 1) = jj^. 

Lattice models : 

(a) (weakly-up/weakly-right) Fix q^, qj > such that qiqj E [0, 1), i = 1, • • • , M, j — 1,- ■ ■ ,N. Take 
X{i,j) ^ giqiQj)- Then 

AI N M N 

^w=nn(i+*^)(i+*/^)' ZM,N^i[i[{i~mjr'. (7.3) 

i—lj — l z— Ij — 1 

(b) (weakly-up/strictly-right) Fix qi,qj > 0, i — 1,- ■ ■ , M, j — I,--- ,iV. Take X{i,j) ^ b{qiqj). 
Then 

M N M N 

^w=nn(i+*^)(i-'?,/^)~'' ZALN^YiY{{i+mj). (7.4) 

(c) (strictly-up/strictly-right) Fix qi, qj > such that q^qj e [0, 1), i = 1, • • • , M , j = 1, • • • ^N. Take 
X{i,j) ^ g{qiqj)- Then 

M N M N 

V{^)=X{X{{^-q^^r\^-q3/^r\ ZALN=X{X{{l~mi)-'- (7.5) 

i—lj—l i—lj—l 

Lattice-Continuum models : 

19 



(d) (Poisson : weakly-right) Fix qi > 0, i = I,--- ,iV. Consider a Poisson process of rate qi on 
M X {i} c M X {1, • • • , N}. Let L{t, N) be the length of the longest weakly-up/weakly-right path 
from (0,1) to {t,N). Then F{L{t,N) < £) = ■^E{jf,u(^e){ip) with 

N N 

^{z) = e''l[{l + qjz), Z,,^ = []e*«'. (7.6) 

1=1 4=1 

(e) (Poisson : strictly-right) In (f), replace the admissible path as weakly- up/strictly-right paths. 
Then 

Af N 

4 = 1 j=l 

Continuum model : 

(f) (Poisson process) This is the model in the square case of Section |[ We have F{L{t) < i) = 
■j^^ueuwif) with 

<^(z)=e*(^+^"\ Zt^e-'\ (7.8) 

Remark 7.1. In (d) and (e), weakly-up condition can be changed to strictly-up, since in a Poisson 
process, the event of two points occurring at the same position has measure 0. 

Remark 7.2. We could have mixture of the above models. Then the corresponding ip is simply multipli- 
cation of each symbol. See |Sq/ for the most general setting. 

7.2 Symmetrized models 

For models (a), (c), (f), one can think of 4 different symmetrized models (see section 7 of p)- We 
could impose that the model is symmetric under (1) rotation about the center, (2) reflection about the 
diagonal, (3) reflection about the anti-diagonal and (4) reflection about both the diagonal and the anti- 
diagonal. In each case, there are still determinant formulas for the distribution of the longest up/right 
path. However, for the symmetry types (2), (3), (4) for (a), (c), we only have formulas for the square 
model M = N . Depending on the symmetry type, instead of Toeplitz determinant above, we sometimes 
have Hankel determinant. 

Here we will consider only the case (2). That is, for example in (a), we impose that X{j,i) should 
be equal to the value oi X(i,j) for each i,j. In each case, the algebraic formula is of the form 

P(i < ^) = I ^ueoiO det(^((7)). (7.9) 

Note that now we have expectation over the orthogonal group. We will specify the function ifj{z). 
Let g'{a, q) be the random variable given by 

'(«,9) = fc)-^^a'^'"°''V. (7.10) 
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(a-S) Fix qi>Q,i = l,---,N and a > such that qi E [0, 1) and aqi E [0, 1). Take X{i, j) = X{j, i) ^ 
giqtqj) for i ^ j and X{i,i) ~ g{aqi). Then the longest weakly-up/weakly-right L{N\a) path 
from (1, 1) to {N,N) satisfies (^ with 



N N 

il^{z)^{l + az)\{{l + q,z), Z = Y{{l-aq,r^ [] {I - q,q,)-\ (7.11) 

1=1 1=1 l<i<j<N 

(c-S) Fix qi > {), i = \,- ■ ■ ,N and a > such that qi E [0,1). Take X{i,i) ~ X{j,i) ~ giqtqj) for 
i ^ j and X(i, i) ^ g'{a, qi). Then the longest strictly-up/strictly-right L{N; a) path from (1, 1) 



to {N,N) satisfies (7.9) with 



N N 

^{z) = {l + az)'[[{l-q,zr\ Z = l[{l + aq,){l-qfr' [] {1 - q,q,)-\ (7.12) 

i=l i=l l<i<j<JV 

(f-S) The Triangle case of Section 0. Take a Poisson process of rate 1 in the region x < y in the 
xy-plane. Then take the mirror image of the points about the line x ^ y. In addition, take a 
Poisson process of rate a > on the line x = y. Then the longest up/right path L(i; a) from 



(0, 0) to (i, t) satisfies ([TJD with 

V'(^) = (l + a^)e*^ Z = e"*+*'/2. (7.13) 
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